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ABSTRACT 

A numerical parametric study is conducted to simulate shock-induced combustion under 
various free-stream conditions and varying blunt body diameter. A steady combustion front is 
established if the free-stream Mach number is above the Chapman-Jouguet speed of the mixture, 
whereas an unsteady reaction front is established if the free-stream Mach number is below or at 
the Chapman-Jouguet speed of the mixture. The above two cases have been simulated for Mach 
5.11 and Mach 6.46 with a projectile diameter of 15 mm. Mach 5.11, which is an underdriven 
case, shows an unsteady reaction front, whereas Mach 6.46, which is an overdriven case, shows a 
steady reaction front. Next for Mach 5.11, reducing the diameter to 2.5 mm causes the instabilities 
to disappear, whereas for Mach 6.46, increasing the diameter of the projectile to 225 mm causes 
the instabilities to reappear, indicating that Chapman-Jouguet speed is not the only deciding 
factor for these instabilities to trigger. The other key parameters are the projectile diameter, 
induction time, activation energy and the heat release. The appearance and disappaerance of the 
instabilities have been explained by the one-dimensional wave interaction model. 
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INTRODUCTION 


Some of the applications of supersonic research are the proposed National Aero-Space Plane, 
the ram accelerator and Trans Atmospheric Vehicles. All of these vehicles rely on air-breathing 
propulsion. The air-breathing engine removes the requirement to carry an oxidizer inside the 
vehicle. Further, there is a substantial savings in weight and, therefore, the payload is higher. An 
efficient propulsion system at hypersonic speeds requires that combustion take place at supersonic 
speeds, i.e., combustion at supersonic mean airstream speed. The supersonic combustion speed is 
required for airbreathers at high Mach numbers to reduce the unacceptable losses associated with 
inlet deceleration to subsonic conditions. The scramjet [1,2] is an integrated airframe-propulsion 
concept for a hypersonic airplane. Significantly increased propulsion efficiency in the Mach 
6 to 15 range may be provided by the conventional diffusion-burning, air-breathing scramjet. 
For airbreathers above Mach 15, another propulsion approach is required, e.g., greatly enhanced 
(reduced loss) conventional scramjet and detonation wave engines. 

Another proposed hypervelocity air-breathing propulsion scheme is the shock-induced com- 
bustion or oblique detonation wave engine (ODWE [3]) where a shock is employed to increase 
the temperature of the premixed fuel and air to a point that creates chemical reaction. The 
term detonation is applied to the process where a shock and reaction front follow each other 
very closely and are pressure coupled, while shock-induced combustion implies that the shock 
wave and reaction front are decoupled. Thus, detonation is a limiting case of shock-induced 
combustion. Up to about Mach 15, a scramjet has better performance than an ODWE, but after 
that, ODWE performs better. This engine requires even more extreme integration with the “air- 
frame.” To avoid preignition, the fuel is injected from the vehicle forebody/nose region, inside 
the shock layer but outside the boundary layer. The advantages of a shock-induced combustion 
propulsion system include lower inlet losses and significant reduction in combustion weight, size, 
heat transfer, and skin friction losses and cooling requirements. These advantages are due to 
forebody region fuel-air mixing which obviates the need to contain the fuel-air mixing process 
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within a long combustor at high pressure. 


Ballistic range experiments performed in the 1960’s and 1970’s provide an excellent source 
for studying supersonic combustion/detonation. The physics of these ballistic range flows are 
predominantly driven by reaction kinetics and convection phenomena. Thus, the complications 
and uncertainties of diffusion and mixing are removed from the problem. Any discrepancy 
between the experimental data and numerical calculations can be attributed to either numerical 
errors or the improperly modelled chemical kinetics. Zeldovich [4] suggested that combustion can 
be stabilized by the shock wave produced by bodies moving at supersonic speeds in combustible 
mixtures at various level of overdrive. In these experiments, projectiles were fired in different 
premixed fuel-air mixtures, and detonation structures around the projectiles were recorded. Every 
gas mixture has a detonation wave velocity known as the Chapman- Jouget (C-J) velocity, which 
is characteristic of the mixture. The detonation wave velocity is known as the C-J velocity of the 
mixture whenever the normal component of the flow velocity following the detonation wave is 
sonic. If the normal component of the flow velocity is subsonic, however, the detonation wave 
velocity is called overdriven and if supersonic, then it is known as underdriven. The free-stream 
velocity is referred to as superdetonative if the free-stream velocity of the projectile is above 
the C-J velocity of the reactive mixture. The detonation wave structure is highly unstable for 
projectile velocities that are less than the C-J velocity of the mixture. The detonation or reaction 
front structure shows a coupled shock-deflagration system near the stagnation line of the body if 
the projectile is flying above the C-J velocity of the gas mixture. These two fronts separate from 
each other as "one moves away from the stagnation line. The separation between the two fronts 
occurs as soon as the velocity component normal to the bow shock is equal to the detonation 
velocity. The separation between the bow shock and the reaction front is called the induction 
zone. Lehr’s [5] experimental ballistic range shadowgraphs for Mach 5.11 and Mach 6.46 are 
shown in Figs. 1 and 2, respectively. A free-stream temperature of 292 K and a pressure of 
42663.2 N/m^ (320 mm of Hg) are used along with a stoichiometric mixture of hydrogen and air 
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in both cases. In Fig. 1 the projectile velocity is Mach 5.1 1 (which is also the C-J velocity of the 
mixture). It shows separation of two discontinuities by induction zone. The outer discontinuity 
is the bow shock, and the inner discontinuity is the reaction front. The separation between the 
two shocks is minimum near the stagnation region and increases downstream. The induction 
distance is increased as the shock wave becomes weak away from the stagnation region and the 
post-shock temperature reduces. Another feature to be noticed is the presence of instabilities 
in the reaction front. The frequencies of these pulsations were determined to be 1.96 MHz. 
It is seen in Fig. 2 that the reaction front for the Mach 6.46 case is coupled with the shock 
near the stagnation line. The coupling continues until about 60° from the stagnation line, where 
the reaction front starts decoupling from the bow shock. This coupling is caused by very high 
post-shock temperatures, that cause a decrease in the induction distance. A close examination 
of the shadowgraphs reveals that there is an increase in density as the flow crosses the bow 
shock. The increase is noted by the color changes from light to dark. As the flow crosses the 
reaction front, however, the color changes from dark to light, indicating a decrease in density 
across the reaction front. This decrease is due to a large release of energy across the reaction 
front that causes an increase in the temperature; since the pressure remains relatively constant, 
the density must decrease. 

By means of schlieren photography, Reugg and Dorsey [6] investigated the problems and 
effects of stabilizing combustion on a 20 mm diameter spherical projectile in a stoichiometric 
mixture of hydrogen and air at rest. Combustion produced detectable effects on the shapes and 
positions of shock waves at Mach numbers between 4 and 6.5 and above a pressure of one-tenth 
atmosphere. Strong combustion-driven oscillations were observed in front of the sphere with 
frequencies of one-tenth megacycle per second. These oscillations were observed when the Mach 
number was less than 6 for the one-half atmosphere pressure, and less than 5 for the one-quarter 
atmosphere pressure. They also reported the results with methane-air and pentane-air mixtures. 

Similar experiments were performed by Behren et al. [7] by firing 9 mm diameter plastic 
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spheres into hydrogen-air and hydrogen-oxygen mixture at velocities of 1500-3000 m/sec. 
Similar behaviour of transition from detonation wave stability to instability was observed when 
the projectile velocities were decreased to nearly C-J velocity of the mixture. 

Chemyi [8] reported on the experiments of firing a sphere-nosed cylinder into hydrogen- 
oxygen and hydrogen-air mixtures. As in the case for self-sustained detonation, the flow field is 
highly unstable under a wide range of conditions. In the case of blunt projectiles and exothermic 
gas mixtures, the instabilities generate a highly regular, periodic flow structure under certain 
conditions. He suggested the similarity of these instabilities to the cellular structure of flame. 

McVey and Toong [9] conducted similar experiments where projectiles were fired into 
lean acetylene-oxygen and stoichiometric hydrogen-air mixtures. They developed the wave 
interaction model to explain the detonation wave structure instabilities. Their model explains 
how compression waves can be formed when a new reaction front develops in the induction zone 
between the normal segment of the bow shock and the original reaction front. These compression 
waves lead to a cyclic process which is compatible with most of the observed features of the 
flow. However, the compression wave strength remained unresolved in their wave-interaction 
model, which is an important factor in determining if such a model is physically possible. Alpert 
and Toong [10] included the effect of the compression wave strength and proposed a modified 
form of the wave-interaction model. 

Several researchers [11-15] have recently attempted to numerically simulate Lehr’s ballistic 
range experiments [5]. Youngster et al. [11] and Lee and Deiwart [12] simulated Lehr’s 
experimental data for Machs 4.18, 5.11, and 6.46. They used Euler equations coupled with 
species equations to capture the shock and the reaction front. The reaction model used was 
a hydrogen-air mixture of six species and an inert gas such as Argon or Nitrogen and eight 
reactions. The flow field was found to be steady despite the experimental evidence that the flow 
field is unsteady. For the test conditions of stoichiometric hydrogen-air mixture, the detonation 
wave speed of the mixture is Mach 5.11. It has been demonstrated experimentally by Lehr that 
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Machs 5.11 and 4.18 show detonation wave structural instabilities that disappear if the flight 
Mach number is increased beyond Mach 5.11. Further, the flow field was not well resolved. 
Their blunt body calculations used 32x32 and 57x41 size grids, respectively. These grids were 
not sufficient to resolve the flow field correctly. 

Wilson and MacCormack [13] conducted a detailed numerical investigation of the shock- 
induced combustion phenomena. Euler equations and a 13-species, 33-reactions chemistry 
model was used. The validity of the reaction models and the importance of grid resolution 
needed to properly model the flow physics were also shown. Highly resolved calculations for 
Lehr’s Mach 5.11 and Mach 6.46 cases with an adaptive grid were performed. The calculations 
were not time accurate; therefore, the unsteady behavior was not captured. 

Sussman and Wilson [14] also studied the instabilities in the reaction front for a Mach 
number of 4.79. Euler equations and a 13-species, 33-reactions chemistry model was again 
used. They have proposed a new formulation based on logarithmic transformation. The number 
of grid points needed to properly resolve the reaction front is greatly reduced. They successfully 
simulated the unsteady case. However, the frequency was slightly underpredicted. 

Matsuo and Fujiwara [15] have studied the instabilities of shock-induced combustion around 
an axisymmetric blunt body with Euler equations and a simplified two-step chemistry model. The 
growth of periodic instabilities by a series of simulations with various tip radii was investigated; 
these periodic instabilities are related to shock-standoff distance and induction length. They 
proposed a new model based on McVey and Toong’s model [9] that explained the instabilities 
in the reaction front. 

The appearance or disappearance of the instabilities in the structure of the reaction front 
with various Mach numbers has been studied in details, and a satisfactory explanation of the 
various flow phenomena have been explained by the wave-interaction model by McVey and 
Toong [9]. The instabilities originate in the induction zone that separates the bow shock and 
the exothermic reaction front in the nose region of the flow field and then spreads outwards. 
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But the disappearance of instabilities by reducing the projectile diameter while keeping the same 
Mach number has not been clarified yet. The objective of this study is to isolate the various 
key parameters for the appearance and disappearance of these instabilities such as projectile 
diameter, induction time, heat release, activation energy and to explain the physics of the 
observed phenomena by a wave interaction model. In order to capture the physical instabilities, 
the calculations must be carried out for long times to ensure that all relevant time scales are being 
captured. Since all numerical schemes have some numerical diffusion, which is dependent on 
the grid resolution, a coarse grid may damp these oscillations. Further, the numerical damping 
added to the scheme in the reaction front vicinity may damp or alter the instability modes. The 
axisymmetric version of the SPARK2D code [16], which incorporates a 9-species, 18-reactions 
combustion model for hydrogen-air mixtures, is used to carry out the analysis. 

BASIC GOVERNING EQUATIONS 

The physical model for analyzing the flow field is described by the Navier-Stokes and species 
continuity equations. For two-dimensional axisymmetric flows, these equations are expressed 
in physical coordinates as [16] 


^ ^ ^ 

dt dx dy 


( 1 ) 


where vectors U, F, G, and H are written as 


U = 


P 

pu 

pv 

pE 

Pfi 


6 




7 



The other terms appearing in vectors F, G, and H are defined as 
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In Eq. (1) only (Ns— 1) species equations need to be considered in the formulation since 
the mass fraction of the species is prescribed by satisfying the constraint equation 

E f 1 

h = 1 (12) 

t=l 


The specific heat at constant pressure for each species is prescribed in Eq. (1 1) by a fourth-order 
polynomial in temperature. The diffusion velocity of the species is obtained by solving the 
Stefan-Maxwell equation, neglecting body force and thermal diffusion effects. 




(14) 


This equation has to be applied only to (Ns— 1) species. The diffusion velocity for the remaining 

. . . . . . ^ 
species IS prescribed by satisfying the constraint equation ^ fiVi = 0, which ensures the 

i=l 

consistency. 


CHEMISTRY AND THERMODYNAMIC MODELS 

Chemical reaction rate expressions are usually determined by summing the contributions 
from each relevant reaction path to obtain the total rate of change of each species. Each path is 
governed by a law of mass action expression in which the rate constants can be determined from 
a temperature dependent Arrhenius expression. In vector H, the term ui — MiCi represents the 
net rate of production of species i in all chemical reactions and is modelled as follows : 


N, N, 




"-S- • 

Jt‘^1 J 

; j = 1, 
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where Eq. (14) is a representation of an N^-step chemical reaction and Eq. (15) is the production 
rate for the species, as determined from the law of mass action. The reaction constants /c^ 
and «bj are calculated from the following equations 
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(18) 


Kf. 

Kbj = 

l^eqj 

The equilibrium constant appearing in Eq. (18) is given by 

where 

N, N, 

i=l i=l 

Ns Ns 

t=l 1=1 

| = A.r(i-/„D-(f)r^-(|)r^ 

The forward rate for each reaction is determined by Eq. (17) which is based on the Arrhenius 
law. The appropriate constants Aj,aj, and ej for the H 2 -air reaction system can be found in 
[17]. The reverse rate is then calculated from Eq. (18). The Gibb’s free energy of each species 
in gi is obtained from the expression for Cpi. 

The hydrogen-air combustion mechanism used in this work is based on the Jachimowski 
hydrogen-air model [17] which uses 9 species and 18 reactions. The species are N 2 , O 2 , H 2 , 
OH, H, O, H 2 O, HO 2 , and H 2 O 2 . Each of the 18 reactions can proceed in the forward and 
backward directions. The reactions are 

1) O 2 + H 2 ^OH + OH 

2) O 2 + H^ OH + O 

3) H 2 + OH ^ H 2 O + H 

4) H 2 + O ^ OH +H 

5) OH + OH ^ H 2 O + O ■ 

6) OH + H + M ^ H 2 O + M 

7) H + H + M^H2+M 
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(19) 

( 20 ) 
( 21 ) 

( 22 ) 



8 ) H + ©2 + M ^ HO2 + M 

9 ) HO2 + OH ^ H2O + M 

10 ) HO2 + M ^ H2 + O2 

11 ) HO2 + H OH + OH 

12 ) HO2 + O ^ OH + O2 

13 ) HO2 + HO2 ^ H2O2 + O2 

14 ) HO2 + H2 ^ H2O2 + H 

15 ) H2O2 + OH ^ H2O + HO2 

16 ) H2O2 + H ^ H2O + OH 

17 ) H2O2 + O ^ HO2 + OH 

18 ) H2O2 + M ^ OH + OH + M 

METHOD OF SOLUTION 

The governing equations are transformed from the physical domain (x, y) to a computational 
domain (^, rj) using an algebraic grid generation technique. In the computational domain, Eq. 
( 1 ) is expressed as 


^ ^ ^ - U 

dt 


( 23 ) 


where 

U Fyrj - Gxr, 


G = Gx^~ Fy^, H^HJ 


( 24 ) 


J — ^^yrf 

MacCormack’s [ 18 ] method is used to solve the governing equations. The scheme is second-order 
accurate in time and space, which results in a spatially and temporally discrete, simultaneous 
system of equations at each grid point. The system of equations is solved subject to initial and 
boundary conditions. At the supersonic inflow boundary, all flow quantities are specified as free- 
stream conditions. At the supersonic outflow boundary, all flow quantities are extrapolated from 
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interior grid points. Although full Navier-Stokes (N-S) equations are used, the slip conditions 
are used to numerically simulate the inviscid flow. A flow tangency or slip boundary condition 
is implied on solid wall. The wall temperature and pressure are extrapolated from interior grid 
points. Initial conditions are obtained by specifying free-stream conditions throughout the flow 
field; the resulting set of equations is marched in time. 

The Lax-Wendroff type schemes are inherently unstable and, hence, higher order numerical 
dissipation terms are often necessary to get a stable solution. For a non-reacting flow field, an 
artificial viscosity based on temperature and/or pressure is traditionally used, but in chemically 
reacting flows, in addition to temperature and pressure gradients, there may be strong species 
concentration gradients. To suppress the numerical oscillations in the induction zone (where 
the gradients in the concentration of reactants and products are very strong), additional artificial 
viscosity based on H 2 O mass fraction similar to the one used by Singh et al. [19] is used. 


RESULTS AND DISCUSSIONS 

The key parameters for the onset of periodic unsteadiness have been identified as (1) 
induction time, (2) reaction rate constant, (3) activation energy, (4) heat release and (5) projectile 
nose radius. In this study we shall be discussing the effect of various nose radii on the stability 
of the reaction front while keeping the first four parameters constant by choosing a particular 
reaction model and by fixing the free-stream Mach number. 

When a blunt body is moving through a reactive mixture at hypersonic speeds, a bow shock 
is formed ahead of the body, and the temperature of the fuel-air mixture, after the bow shock, 
is sufficiently high to initiate the reaction. Once ignition starts, chemical energy is released and 
another discontinuity known as the reaction front is formed. In the induction zone, temperature 
and pressure remain relatively constant at the post shock conditions, while the concentrations of 
radicals build up very rapidly. The flow attains equilibrium due to large residence time in the 
stagnation zone, while away from the stagnation zone, the flow is in a state of non-equilibrium. 
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The numerical simulation is carried out for the free-stream conditions given in Table 1 

Table 1 

Free-stream conditions and different nose radii used 


IVIoo 

Dp, mm 

Poo, N/m^ 

Too, k 

5.11 

i) 15 

42663.22 (320 mm) 

292 


ii) 2.5 



6.46 

iii) 15 

42663.22 (320 mm) 

292 


iv) 150 




v) 225 




The premixed fuel oxidizer mixture is taken as 2 H 2 + O 2 + 3.76N2, and the stoichiometric 
chemical reaction for the system can be written as 

2H2 + O2 + 3.76N2 ^ 2H2O + 3.76N2 

Calculations have been carried out for a grid with 197 points in the circumferential direction 
and 152 points in the normal direction. This grid was chosen based on the earlier work by Ahuja 
et al. [20] where the flow field was shown to be adequately resolved with this grid. For the 
present stoichiometric hydrogen-air mixture, the C-J velocity is Mach 5.1 1. Thus, for the Mach 
6.46 case, the projectile speed is significantly above the detonation velocity of the mixture. If 
the free-stream velocity of the projectile is around the C-J detonation velocity of the mixture, 
unsteady flow phenomenon can occur. As given in Table 1, we shall be discussing five cases, 
two cases with Mach 5.11 and projectile diameters of 15 mm and 2.5 mm, and the remaining 
three cases with a Mach number of 6.46 and projectile diameters of 15 mm, 150 mm, and 225 
mm. In all the calculations the residuals dropped by three orders in 12,000 iterations and then 
remained constant. 
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Case (i): Mach 5.11 and projectile diameter of 15 mm 

First the projectile diameter was chosen as 15 mm for the free-stream Mach number of 
5.11. Figure 3 shows the contour plot of temperature, and Fig. 4 shows the contour plot of 
water mass fraction. The bow shock is very smooth, but the reaction front shows oscillations, 
which is more clearly seen in the inset enlarged view. The bow shock and the reaction front 
are separated from each other by the induction distance as is clearly seen in Fig. 5 which is the 
contour plot of density. The separation is minimum at the stagnation line and increases away 
from it because of low post-shock temperature away from the stagnation region. A comparison 
with Fig. 1 shows that all flow features are very well captured. The calculated frequencies of 
these oscillations [20] were found to be in good agreement with the experimentally observed 
frequency [5, 21]. 

Case (ii): Mach 5.11 and projectile diameter of 2.5 mm 

Next the diameter of the projectile was reduced to 2.5 mm while keeping the same free-stream 
Mach number of 5.1 1. Other free-stream conditions were also kept the same. Figure 6 shows the 
temperature contours, and Fig. 7 shows the density contours. Both the shock and the reaction 
front are very smooth. This is much more clear from the water mass fraction contours shown in 
Fig. 8 with an inset enlarged view showing smooth reaction front. Thus, reducing the projectile 
diameter caused the instabilities to disappear. 

Case (Hi): Mach 6.46 and projectile diameter of 15 mm 

In this case the incoming Mach number was increased to 6.46, making it a superdetonative 
case, while keeping the projectile diameter of 15 mm as in case (i). Figure 9 shows the contour 
plot of density for Mach 6.46. The bow shock and the reaction front can be clearly seen in the 
figure. They are coupled with each other near the stagnation line and up to about 60 degrees 
from the nose, at which point they start decoupling from each other by the induction distance. 
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This occurs because bow shock is almost normal near the stagnation line and the post-shock 
temperature is maximum. For Mach 6.46, a very small induction distance occurs as a result of 
the post-shock temperature remaining significantly high up to some distance near the stagnation 
zone. Away from the stagnation line, the induction distance is increased as a result of decreasing 
shock strength and post-shock temperature. A comparison with Fig. 2 shows that all the flow 
features are very well captured. Figure 10 shows the contour plot of ternperamre. A comparison 
with Fig. 3 shows that when the Mach number is increased from a C-J Mach number of 
5.11 to a superdetonative Mach number of 6.46, while keeping the same projectile diameter, 
the instabilities in the reaction front disappear. Thus, for the Mach 6.46 case with projectile 
diameter of 15 mm, both the bow shock and the reaction front have a smooth profile. 

Case (iv): Mach 6.46 and projectile diameter of 150 mm 

Next the projectile diameter was increased to 150 mm while keeping the same Mach number 
of 6.46. Figure 1 1 shows the enlarged view of contour plot of water mass fraction for Mach 
6.46, but with projectile diameter of 150 mm. The figure shows the waviness of the reaction 
front, indicating that the reaction front is turning towards instability. 

Case (v): Mach 6.46 and projectile diameter of 225 mm 

Next, the diameter of the projectile was increased to 225 mm while keeping the same free- 
stream Mach number of 6.46 as in the previous case. Figure 12 shows the enlarged view of 
the water mass fraction contours, and the periodic instabilities of the reaction front are clearly 
evident. Thus, increasing the projectile diameter from 15 mm to 150 mm and then to 225 mm, 
while keeping the same superdetonative Mach number of 6.46, causes the transition of a stable 
reaction front to an unstable one. .The results of the above five cases have been summarized 
in Table 2. 
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Table 2 


Summary of five test cases 


Moo 

Projectile Diameter 

Profile of Reaction Front 


D 


5.11 

(i) 15 mm 

Unstable 


(ii) 2.5 mm 

Stable 

6.46 

(iii) 15 mm 

Stable 


(iv) 150 mm 

Slightly Unstable 


(v) 225 mm 

Unstable 


In order to have a clear understanding of the origin and propagation of the instabilities, an 
x-t diagram for water mass fraction on the stagnation streamline between the bow shock and 
the reaction front is drawn. 

Figure 13 shows the x-t plot of water mass fraction along the stagnation line for the Mach 
5.11 for a 15 mm projectile diameter with an overlay of pressure to show the location of shock 
front. Fig. 14 shows the x-t plot of water mass fraction along the stagnation line for Mach 5.1 1, 
but with a projectile diameter of 2.5 mm. In the former case the reaction front clearly shows 
periodic oscillations, whereas the later case shows a smooth reaction front. The instabilities for 
Mach numbers lower than the C-J Mach number are due to the ignition delay. What causes these 
instabilities to disappear for the same Mach number of 5.1 1 but lower projectile diameter will be 
explained in subsequent paragraphs. A similar trend which was observed for the superdetonative 
case of Mach 6.46, where increasing the diameter causes the instabilities to reappear, shall also 
be discussed. A qualitative study of the frequencies, shock structure, pressure, and stagnation 
temperature associated with the projectile velocities, shows that the reciprocal frequencies or 
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periods of oscillation are equal to the induction time for the ignition of hydrogen-air mixtures. 

Many features of the non-steady flow regimes appear to originate from near the stagnation 
region behind the normal segment of the bow shock. Also, there is clearly a relationship 
between the periodicity of these flow features and the magnitude of the chemical induction time. 
The appearance and disappearance of these instabilities with various nose diameters have been 
explained with a wave interaction model as proposed by McVey and Toong [9]. Figure 15 
shows a schematic x-t diagram for the wave interaction model for Mach 5.1 1 and with projectile 
diameters of 15 mm and 2.5 mm. For the diameter of 15 mm, the beginning of the cycle is shown 
at a time when the overtaking of the bow shock by a compression wave in the induction zone 
takes place. Since chemical induction time is exponentially dependent on the gas temperature, 
the contact discontinuity (DCi) generated when the incident compression wave overtakes the 
bow shock has an important effect on the location of a reaction front. Exothermic reaction 
takes place on the upstream side of the contact discontinuity (DCi) before beginning on the 
downstream side due to the difference in the gas temperature on the upstream and downstream 
sides of the contact discontinuity. It is this beginning of exothermic reaction on only one side 
of the contact discontinuity which constitutes the formation of a new reaction front. The ratio 
of chemical induction time across the discontinuity is a measure of where the new reaction front 
will be located relative to the original reaction front. Since the gases are hot on the upstream side 
of the contact discontinuity, they begin to react, generating compression waves which propagate 
upstream and downstream. This is shown in the figure by (CW 2 ) and (CW 3 ). At a somewhat later 
time, the contact discontinuity reaches the position of the original reaction front, extinguishing 
the reaction at this front and generating rarefaction waves (RWi) and (RW 2 ). The upstream- 
facing compression wave (CW 2 ) eventually overtakes the bow shock and thus results in the 
generation of another contact discontinuity (DC 2 ). Presence of both rarefaction and compression 
waves in the induction zone makes it possible for a periodic wave interaction process. 

It is clear from the above discussion that a compression wave moving towards the bow shock 
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is overtaken by a rarefaction wave in each cycle of the wave-interaction process. During this 
process the compression wave gets attenuated by a rarefaction wave. If Msi is the Mach number 
of the incident compression wave (CWi) or weak shock, then (MsiW is generally of the order 
of 1.01. A wave interaction process cannot be initiated by an unattenuated compression wave 
since the Mach number ( Mgi) of such a wave is much greater than (Msi)max- On the other hand, 
if a sufficiently weak compression wave does initiate the wave interaction process, the attenuation 
of the first and each subsequent compression wave by a rarefaction wave must be sustained in 
order that the Mach number of each reaction shock, which finally overtakes the bow shock, be 
less than (Msi)max. When this condition is satisfied, then only the cyclic process is feasible. The 
validation of this phenomena with a shock-fitting method shall be presented in future studies. 

The importance of a compression wave reflecting off the projectile nose has not been clarified 
in the McVey and Toong model. The compression wave which travels towards the projectile 
reflects from it and finally overtakes the bow shock. This overtaking of the bow shock by 
the reflected compression wave from the projectile takes place almost at the same time that 
a newly formed compression wave from the new reaction front has overtaken the bow shock. 
Consequently, this contributes to strengthening the compression wave coming from the reaction 
front. 

Thus, when we compare the unstable Mach 5.11 case with a projectile diameter of 15 mm 
with the stable case of a projectile diameter of 2.5 mm, it is clear that the 15 mm diameter case 
has greater shock stand-off distance than that with a diameter of 2.5 mm. Thus, the reflected 
compression wave from the projectile nose coming back to overtake the bow shock to strengthen 
the compression wave originating at the reaction front, would have become quite weak. This 
is because it undergoes a large number of interactions with other compression waves between 
the reaction front and the projectile nose. As shown in the figure, the strength of the shock 
P 5 /P 3 is less than the strength of the original reflected shock which is P 2 /P 1 . Now, if the shock 
standoff distance is larger, the reflected compression wave would encounter a greater number 


18 



of interactions, and with every interaction, its strength would decrease. Although this reflected 
compression wave strengthens the compression wave originating at the reaction front, the Ms 
is still less than (Msi)max &nd, therefore, the periodic instabilities are sustained. On the other 
hand, when the projectile diameter is 2.5 mm, the reflected compression wave has to travel 
a very short distance before it reaches the reaction front to strengthen the compression wave 
originating there. Moreover, because of the short shock stand-off distance, it does not encounter 
many interactions with other compression waves. Although P 3 /P 2 of the reflected wave is less 
than P 2 /P 1 of the compression wave reaching the projectile, it is still very strong by the time it 
reaches the reaction front. Thus, when it overtakes the bow shock, it strengthens the compression 
wave much more than the preceeding 15 mm diameter case. This causes Mg of the compression 
wave to be much greater than (Msi)max> and this leads to unattenuated compression waves and, 
hence, periodicity disappears. 

Again, when referring to the Mach 6.46 case with a projectile diameter of 15 mm, a stable 
reaction front is observed, whereas increasing the diameter to 225 mm causes the instabilities 
to reappear. For the 225 mm case where the shock stand-off distance is much higher than the 
15 mm case, the shock wave reflected off the projectile body becomes much weaker before 
strengthening the compression wave generated at the new reaction front. Thus, the compression 
wave generated at the new reaction front has Mg less than (Mgi)max and, therefore, a cyclic 
process is sustained. On the other hand, for the projectile diameter of 15 mm where the shock 
stand-off distance is much smaller, the reflected compression wave from the projectile body 
strengthens the compression wave generated at the new reaction front much more and, therefore, 
leads to unattenuated compression waves. Consequently, the cyclic process for the 15 mm 
diameter and Mach 6.46 case cannot be sustained. 
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CONCLUSIONS 


Numerical studies with different projectile diameters and Mach numbers indicate that 
Chapman-Jouget velocity is not the only parameter for triggering the instabilities. The other 
parameters are induction time, reaction rate constant, projectile diameter, activation energy and 
heat release. The reason for the appearance and disappearance of instabilities could be attributed 
to attenuated or unattenuated compression waves which depend upon the strength of the reflected 
compression wave from projectile body, the strength of which in turn depends upon the shock 
stand-off distance and, hence, the projectile diameter. Thus, it is concluded that a subdetonative 
case can also lead to a stable reaction front by having an appropriate small diameter projec- 
tile, and a superdetonative case can also lead to an unstable reaction front by having a larger 
diameter projectile. 
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Figure 1 Shadowgraph of a spherical nose projectile moving at 
Mach 5.11 into a premixed stoichiometric hydrogen-air mixture. 



Figure 2 Shadowgraph of a spherical nose projectile moving at 
Mach 6.46 into a premixed stoichiometric hydrogen-air mixture. 
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Figure 3 Contour plot of temperature for Mach 5.11 and projectile diameter 15 mm. 
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Figure 4 Contour plot of water mass fraction for Mach 5.11 and projectile diameter 15 mm. 
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Figure 5 Contour plot of density for Mach 5.11 and projectile diameter 15 mm. 
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Figure 6 Contour plot of temperature for Mach 5.11 and projectile diameter 2.5 mm. 
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Figure 7 Contour plot of density for Mach 5.11 and projectile diameter 2.5 mm. 
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Figure 8 Contour plot of water mass fraction for Mach 5.11 and projectile diameter 2.5 mm. 
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Figure 9 Contour plot of density for Mach 6.46 and projectile diameter 15 mm. 
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Figure 10 Contour plot of temperature for Mach 6.46 and projectile diameter 15 mm. 
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Figure 1 1 Enlarged view of contour plot of water mass 
fraction for Mach 6.46 and projectile diameter 150 mm. 
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Figure 12 Enlarged view of contour plot of water mass 
fraction for Mach 6.46 and projectile diameter 225 mm. 
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Figure 13 x-t Plot of water mass fraction along stagnation 
streamline for Mach 5.11 and Projectile diameter 15mm 
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Figure 14 x-t Plot of water mass fraction along stagnation 
streamline for Mach 5.11 and Projectile diameter 2.5 mm 
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ABSTRACT 

A numerical parametric study is conducted to sim- 
ulate shock-induced combustion under various free- 
stream conditions and varying blunt body diameter. A 
steady combustion front is established if the free-stream 
Mach number is above the Chapman-Jouguet speed of 
the mixture, whereas an unsteady reaction front is es- 
tablished if the free-stream Mach number is below or at 
the Chapman-Jouguet speed of the mixture. The above 
two cases have been simulated for Mach 5.1 1 and Mach 
6.46 with a projectile diameter of 15 mm. Mach 5.11, 
which is an underdriven case, shows an unsteady reac- 
tion front, whereas Mach 6.46, which is an overdriven 
case, shows a steady reaction front. Next for Mach 
5.11, reducing the diameter to 2.5 mm causes the insta- 
bilities to disappear, whereas for Mach 6.46, increas- 
ing the diameter of the projectile to 225 mm causes 
the instabilities to reappear, indicating that Chapman- 
Jouguet speed is not the only deciding factor for these 
instabilities to trigger. The other key parameters are 
the projectile diameter, induction time, activation en- 
ergy and the heat release. The appearance and disap- 
paerance of the instabilities have been explained by the 
one-dimensional wave interaction model. 

NOMENCLATURE 

Aj reaction rate constant for the j* reaction 
Ci concentration of i* species 

Cpi constant pressure specific heat of i‘*’ 

species 

Dij binary diffusion coefficient of i* and j* 
species 

Dp projectile diameter 

E total (internal and kinetic) energy 

€j activation energy of j^ reaction 

hj^ base enthalpy of i^ species 

k thermal conductivity 

Kg forward rate constant for j^’’ reaction 

Kbj backward rate constant for j^** reaction 
Keqj equilibrium constant for j*^ reaction 

Mj molecular weight of i* species 

Ns number of chemical species . 

Nr number of chemical reactions 

p pressure 

Graduate Research Assistant, AIAA Student Member 
^ Eminent Professor, AIAA Associate Fellow 


Ri gas constant of i* species 

T temperature 
u x-component of the velocity 

Ui x-component of the diffusion velocity of 

the i'^ component 

V y'Component of the velocity 

Vi y-component of diffusion velocity of 

i* species 

Xi mole fraction of i* species 
X streamwise coordinate in the physical 
domain 

y normal coordinate in the physical domain 
CTj temperature coefficient in reaction rate 
expression for j'*’ reaction 
AGrj Gibbs free energy change for the j*^ 
reaction 

Anj molar change for the j* reaction 
77 normal coordinate in the computational 

domain 

A second viscosity coefficient 

fi dynamic viscosity 

i/j,- stoichiometric coefficient of reactant 

corresponding to i**’ species and j*** reaction 
Uji stoichiometric coefficient of product 

corresponding to i^^ species and j* reaction 
^ streamwise coordinate in the computational 

domain 
p density 

< 7 x normal stress in the x-direction 

< 7 y normal stress in the y-direction 

Txy shear stress in the xy plane 

a), production rate of i^*’ species 

INTRODUCTION 

Some of the applications of supersonic research 
are the proposed National Aero-Space Plane, the ram 
accelerator and Trans Atmospheric Vehicles. All of 
these vehicles rely on air-breathing propulsion. The 
air-breathing engine removes the requirement to carry 
an oxidizer inside the vehicle. Further, there is a sub- 
stantial savings in weight and, therefore, the payload 
is higher. An efficient propulsion system at hyper- 
sonic speeds requires that combustion take place at su- 
personic speeds, i.e., combustion at supersonic mean 
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airstream speed. The supersonic combustion speed is 
required for airbreathers at high Mach numbers to re- 
duce the unacceptable losses associated with inlet de- 
celeration to subsonic conditions. The scramjet [1,2] 
is an integrated airframe-propulsion concept for a hy- 
personic airplane. Significantly increased propulsion 
efficiency in the Mach 6 to 15 range may be provided 
by the conventional diffusion-burning, air-breathing 
scramjet. For airbreathers above Mach 15, another 
propulsion approach is required, e.g., greatly enhanced 
(reduced loss) conventional scramjet and detonation 
wave engines. 

Another proposed hypervelocity air-breathing 
propulsion scheme is the shock-induced combustion 
or oblique detonation wave engine (ODWE [3]) where 
a shock is employed to increase the temperature of the 
premixed fuel and air to a point that creates chem- 
ical reaction. The term detonation is applied to the 
process where a shock and reaction front follow each 
other very closely and are pressure coupled, while 
shock-induced combustion implies that the shock wave 
and reaction front are decoupled. Thus, detonation 
is a limiting case of shock-induced combustion. Up 
to about Mach 15, a scramjet has better performance 
than an ODWE, but after that, ODWE performs better. 
This engine requires even more extreme integration 
with the “airframe.” To avoid preignition, the fuel is 
injected from the vehicle forebody/nose region, inside 
the shock layer but outside the boundary layer. The 
advantages of a shock-induced combustion propulsion 
system include lower inlet losses and significant re- 
duction in combustion weight, size, heat transfer, and 
skin friction losses and cooling requirements. These 
advantages are due to forebody region fuel-air mixing 
which obviates the need to contain the fuel-air mixing 
process within a long combustor at high pressure. 

Ballistic range experiments performed in the 
1960’s and 1970’s provide an excellent source for 
studying supersonic combustion/detonation. The 
physics of these ballistic range flows are predominantly 
driven by reaction kinetics and convection phenomena. 
Thus, the complications and uncertainties of diffusion 
and mixing are removed from the problem. Any dis- 
crepancy between the experimental data and numerical 
calculations can be attributed to either numerical errors 
or the improperly modelled chemical kinetics. Zel- 
dovich [4] suggested that combustion can be stabilized 
by the shock wave produced by bodies moving at super- 
sonic speeds in combustible mixtures at various level 
of overdrive. In these experiments, projectiles were 
fired in different premixed fuel-air mixtures, and deto- 
nation structures around the projectiles were recorded. 
Every gas mixture has a detonation wave velocity 
known as the Chapman-Jouget (C-J) velocity, which 
is characteristic of the mixture. The detonation wave 
velocity is known as the C-J velocity of the mixture 
whenever the normal component of the flow velocity 


following the detonation wave is sonic. If the normal 
component of the flow velocity is subsonic, however, 
the detonation wave velocity is called overdriven and if 
supersonic, then it is known as underdriven. The free- 
stream velocity is referred to as superdetonative if the 
free-stream velocity of the projectile is above the C-J 
velocity of the reactive mixture. The detonation wave 
structure is highly unstable for projectile velocities that 
are less than the C-J velocity of the mixture. The 
detonation or reaction front structure shows a coupled 
shock-deflagration system near the stagnation line of 
the body if the projectile is flying above the C-J veloc- 
ity of the gas mixture. These two fronts separate from 
each other as one moves away from the stagnation line. 
The separation between the two fronts occurs as soon 
as the velocity component normal to the bow shock 
is equal to the detonation velocity. The separation 
between the bow shock and the reaction front is called 
the induction zone. Lehr’s [5] experimental ballistic 
range shadowgraphs for Mach 5.1 1 and Mach 6.46 are 
shown in Figs. 1 and 2, respectively. A free-stream 
temperature of 292 K and a pressure of 42663.2 N/m^ 
(320 mm of Hg) are used along with a stoichiometric 
mixture of hydrogen and air in both cases. In Fig. 

1 the projectile velocity is Mach 5.11 (which is also 
the C-J velocity of the mixture). It shows separation 
of two discontinuities by induction zone. The outer 
discontinuity is the bow shock, and the inner discon- 
tinuity is the reaction front. The separation between 
the two shocks is minimum near the stagnation region 
and increases downstream. The induction distance is 
increased as the shock wave becomes weak away from 
the stagnation region and the post-shock temperature 
reduces. Another feature to be noticed is the presence 
of instabilities in the reaction front. The frequencies of 
these pulsations were determined to be 1 .96 MHz. It is 
seen in Fig. 2 that the reaction front for the Mach 6.46 
case is coupled with the shock near the stagnation line. 
The coupling continues until about 60® from the stag- 
nation line, where the reaction front starts decoupling 
from the bow shock. This coupling is caused by very 
high post-shock temperatures, that cause a decrease 
in the induction distance. A close examination of the 
shadowgraphs reveals that there is an increase in den- 
sity as the flow crosses the bow shock. The increase 
is noted by the color changes from light to dark. As 
the flow crosses the reaction front, however, the color 
changes from dark to light, indicating a decrease in 
density across the reaction front. This decrease is due 
to a large release of energy across the reaction front 
that causes an increase in the temperature; since the 
pressure remains relatively constant, the density must 
decrease. 

By means of schlieren photography, Reugg and 
Dorsey [6] investigated the problems and effects of 
stabilizing combustion on a 20 mm diameter spheri- 
cal projectile in a stoichiometric mixture of hydrogen 
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and air at rest. Combustion produced detectable effects 
on the shapes and positions of shock waves at Mach 
numbers between 4 and 6.5 and above a pressure of 
one-tenth atmosphere, Strong combustion-driven os- 
cillations were observed in front of the sphere with 
frequencies of one-tenth megacycle per second. These 
oscillations were observed when the Mach number was 
less than 6 for the one-half atmosphere pressure, and 
less than 5 for the one-quarter atmosphere pressure. 
They also reported the results with methane-air and 
pentane-air mixtures. 

Similar experiments were performed by Behren et 
al. [7] by firing 9 mm diameter plastic spheres into 
hydrogen-air and hydrogen-oxygen mixture at veloci- 
ties of 1500—3000 m/sec. Similar behaviour of transi- 
tion from detonation wave stability to instability was 
observed when the projectile velocities were decreased 
to nearly C-J velocity of the mixture. 

Chemyi [8] reported on the experiments of fir- 
ing a sphere-nosed cylinder into hydrogen-oxygen and 
hydrogen-air mixtures. As in the case for self-sustained 
detonation, the flow field is highly unstable under a 
wide range of conditions. In the case of blunt pro- 
jectiles and exothermic gas mixtures, the instabilities 
generate a highly regular, periodic flow structure under 
certain conditions. He suggested the similarity of these 
instabilities to the cellular structure of flame. 

McVey and Toong [9] conducted similar experi- 
ments where projectiles were fired into lean acetylene- 
oxygen and stoichiometric hydrogen-air mixtures. 
They developed the wave interaction model to explain 
the detonation wave structure instabilities. Their model 
explains how compression waves can be formed when 
a new reaction front develops in the induction zone 
between the normal segment of the bow shock and the 
original reaction front. These compression waves lead 
to a cyclic process which is compatible with most of 
the observed features of the flow. However, the com- 
pression wave strength remained unresolved in their 
wave-interaction model, which is an important factor 
in determining if such a model is physically possible. 
Alpert and Toong [10] included the effect of the com- 
pression wave strength and proposed a modified form 
of the wave-interaction model. 

Several researchers [11-15] have recently at- 
tempted to numerically simulate Lehr’s ballistic range 
experiments [5]. Youngster et al. [11] and Lee and 
Deiwart [12] simulated Lehr’s experimental data for 
Machs 4.18, 5.11, and 6.46. They used Euler equa- 
tions coupled with species equations to capture the 
shock and the reaction front. The reaction model used 
was a hydrogen-air mixture of six species and an inert 
gas such as Argon or Nitrogen and eight reactions. 
The flow field was found to be steady despite the 
experimental evidence that the flow field is unsteady. 
For the test conditions of stoichiometric hydrogen-air 


mixture, the detonation wave speed of the mixture is 
Mach 5.11. It has been demonstrated experimentally 
by Lehr that Machs 5.11 and 4.18 show detonation 
wave structural instabilities that disappear if the flight 
Mach number is increased beyond Mach 5,11, Fur- 
ther, the flow field was not well resolved. Their blunt 
body calculations used 32x32 and 57x41 size grids, 
respectively. These grids were not sufficient to resolve 
the flow field correctly. 

Wilson and MacCormack [13] conducted a de- 
tailed numerical investigation of the shock-induced 
combustion phenomena. Euler equations and a 
13-species, 33-reactions chemistry model was used. 
The validity of the reaction models and the importance 
of grid resolution needed to properly model the flow 
physics were also shown. Highly resolved calculations 
for Lehr’s Mach 5.11 and Mach 6.46 cases with an 
adaptive grid were performed. The calculations were 
not time accurate; therefore, the unsteady behavior 
was not captured. 

Sussman and Wilson [14] also studied the insta- 
bilities in the reaction front for a Mach number of 4.79. 
Euler equations and a 13-species, 33-reactions chem- 
istry model was again used. They have proposed a 
new formulation based on logarithmic transformation. 
The number of grid points needed to properly resolve 
the reaction front is greatly reduced. They successfully 
simulated the unsteady case. However, the frequency 
was slightly underpredicted. 

Matsuo and Fujiwara [15] have studied the in- 
stabilities of shock-induced combustion around an ax- 
isymmetric blunt body with Euler equations and a sim- 
plified two-step chemistry model. The growth of peri- 
odic instabilities by a series of simulations with vari- 
ous tip radii was investigated; these periodic instabili- 
ties are related to shock-standoff distance and induction 
length. They proposed a new model based on McVey 
and Toong ’s model [9] that explained the instabilities 
in the reaction front. 

The appearance or disappearance of the instabili- 
ties in the structure of the reaction front with various 
Mach numbers has been studied in details, and a sat- 
isfactory explanation of the various flow phenomena 
have been explained by the wave-interaction model by 
McVey and Toong [9]. The instabilities originate in 
the induction zone that separates the bow shock and 
the exothermic reaction front in the nose region of the 
flow field and then spreads outwards. But the disap- 
pearance of instabilities by reducing the projectile di- 
ameter while keeping the same Mach number has not 
been clarified yet. The objective of this study is to iso- 
late the various key parameters for the appearance and 
disappearance of these instabilities such as projectile di- 
ameter, induction time, heat release, activation energy 
and to explain the physics of the observed phenomena 
by a wave interaction model. In order to capture the 
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physical instabilities, the calculations must be carried 
out for long times to ensure that all relevant time scales 
are being captured. Since all numerical schemes have 
some numerical diffusion, which is dependent on the 
grid resolution, a coarse grid may damp these oscil- 
lations. Further, the numerical damping added to the 
scheme in the reaction front vicinity may damp or alter 
the instability modes. The axisymmetric version of the 
SPARK2D code [16], which incorporates a 9-species, 
18-reactions combustion model for hydrogen-air mix- 
tures, is used to carry out the analysis. 

BASIC GOVERNING EQUATIONS 

The physical model for analyzing the flow field is 
described by the Navier-Stokes and species continuity 
equations. For two-dimensional axisymmetric flows, 
these equations are expressed in physical coordinates 
as [16] 
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dt dx dy 
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where vectors U, F, G, and H are written as 
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In Eq. (1) only (Ns — 1) species equations need to 
be considered in the formulation since the mass fraction 
of the species is prescribed by satisfying the constraint 
equation 

N. 

L/i = l (12) 

i=X 

The specific heat at constant pressure for each species is 
prescribed in Eq. (11) by a fourth-order polynomial in 
temperature. The diffusion velocity of the species 
is obtained by solving the Stefan-Maxwell equation, 
neglecting body force and thermal diffusion effects. 



( 14 ) 
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This equation has to be applied only to (Ns — 1) species. 
The diffusion velocity for the remaining species is pre- 

scribed by satisfying the constraint equation Y!, fi^i = 
0, which ensures the consistency. 

CHEMISTRY AND 
THERMODYNAMIC MODELS 

Chemical reaction rate expressions are usually de- 
termined by summing the contributions from each rel- 
evant reaction path to obtain the total rate of change 
of each species. Each path is governed by a law of 
mass action expression in which the rate constants can 
be determined from a temperature dependent Arrhe- 
nius expression. In vector H, the term uu — MiCi 
represents the net rate of production of species i in all 
chemical reactions and is modelled as follows ; 

N, N. 

= (15) 

i=l t=l 


Nr 


«. = Mi 

j=l 


“ji n - '“i n 


m=l 


m=:l 


(16) 

where Eq. (14) is a representation of an N^-step chem- 
ical reaction and Eq. (15) is the production rate for the 
species, as determined from the law of mass action. 
The reaction constants Kg and /Cbj are calculated from 
the following equations 


Kfj = AjT°‘^ exp 


Rut) ^ 


Kf. 

Kbj — ; j — 1, ....Mr 

The equilibrium constant appearing in Eq. 
given by 
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(18) 
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The forward rate for each reaction is determined 
by Eq. (17) which is based on the Arrhenius law. 


The appropriate constants Aj,aj, andej for the H 2 - 
air reaction system can be found in [17]. The reverse 
rate is then calculated from Eq. (18), The Gibb’s 
free energy of each species in g,- is obtained from the 
expression for Cp,-. 

The hydrogen-air combustion mechanism used in 
this work is based on the Jachimowski hydrogen-air 
model [17] which uses 9 species and 18 reactions. The 
species are N 2 , O 2 , H 2 , OH, H, O, H 2 O, HO 2 , and 
H 2 O 2 . Each of the 18 reactions can proceed in the 
forward and backward directions. The reactions are 

1) O 2 + H 2 ^OH + OH 

2) O 2 + H^ OH + O 

3) H 2 + OH H 2 O + H 

4) H 2 + O ^ OH +K 

5) OH + OH H 2 O + O 

6) OH + H + M H 2 O -t- M 

7) H + H + M^H2 + M 

8) H + O 2 + M ^ HO 2 + M 

9) HO 2 + OH^ H 2 O -H M 

10 ) HO2 + M ^ H2 + O2 

11) HO 2 + H ^ OH + OH 

12) HO 2 + O ^ OH + O 2 

13 ) HO2 + HO2 ^ H2O2 + O2 

14) HO 2 + H 2 ^ H 2 O 2 + H 

15) H 2 O 2 + OH ^ H 2 O -H HO 2 

16) H 2 O 2 + H ^ H 2 O + OH 

17) H 2 O 2 + O ^ HO 2 + OH 

18) H 2 O 2 + M ^ OH + OH + M 

METHOD OF SOLUTION 

The governing equations are transformed from the 
physical domain (x, y) to a computational domain (^, 
77 ) using an algebraic grid generation technique. In the 
computational domain, Eq. (1) is expressed as 


dt drj ~ ^ 

where 

U = UJ,F = Fyr, - Gx, 
G = Gx^~ Fy^,H HJ 
J = 


(23) 


(24) 


MacCormack’s [18] method is used to solve the gov- 
erning equations. The scheme is second-order accu- 
rate in time and space, which results in a spatially and 
temporally discrete, simultaneous system of equations 
at each grid point. The system of equations is solved 
subject to initial and boundary conditions. At the super- 
sonic inflow boundary, all flow quantities are specified 
as free-stream conditions. At the supersonic outflow 
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boundary, all flow quantities are extrapolated from in- 
terior grid points. Although full Navier-Stokes (N-S) 
equations are used, the slip conditions are used to nu- 
merically simulate the inviscid flow. A flow tangency 
or slip boundary condition is implied on solid wall. The 
wall temperature and pressure are extrapolated from 
interior grid points. Initial conditions are obtained by 
specifying free-stream conditions throughout the flow 
field; the resulting set of equations is marched in time. 

The Lax-Wendroff type schemes are inherently 
unstable and, hence, higher order numerical dissipation 
terms are often necessary to get a stable solution. For a 
non-reacting flow field, an artificial viscosity based on 
temperature and/or pressure is traditionally used, but in 
chemically reacting flows, in addition to temperature 
and pressure gradients, there may be strong species 
concentration gradients. To suppress the numerical 
oscillations in the induction zone (where the gradients 
in the concentration of reactants and products are very 
strong), additional artificial viscosity based on H 2 O 
mass fraction similar to the one used by Singh et al. 
[19] is used. 

RESULTS AND DISCUSSIONS 

The key parameters for the onset of periodic un- 
steadiness have been identified as (1) induction time, 
(2) reaction rate constant, (3) activation energy, (4) heat 
release and (5) projectile nose radius. In this study we 
shall be discussing the effect of various nose radii on 
the stability of the reaction front while keeping the first 
four parameters constant by choosing a particular reac- 
tion model and by fixing the free-stream Mach number. 

When a blunt body is moving through a reactive 
mixture at hypersonic speeds, a bow shock is formed 
ahead of the body, and the temperature of the fuel- 
air mixture, after the bow shock, is sufficiently high 
to initiate the reaction. Once ignition starts, chemical 
energy is released and another discontinuity known as 
the reaction front is formed. In the induction zone, 
temperature and pressure remain relatively constant at 
the post shock conditions, while the concentrations of 
radicals build up very rapidly. The flow attains equi- 
librium due to large residence time in the stagnation 
zone, while away from the stagnation zone, the flow is 
in a state of non-equilibrium. The numerical simula- 
tion is carried out for the free-stream conditions given 
in Table 1. 


Table 1 

Free-stream conditions and 
different nose radii used 


Moo 

Dp, mm 

Poo, N/m^ 

Too, k 

5.11 

i) 15 

ii) 2.5 

42663.22 
(320 mm) 

292 

6.46 

iii) 15 

iv) 150 
V) 225 

42663.22 
(320 mm) 

292 


The premixed fuel oxidizer mixture is taken as 
2 H 2 + O 2 + 3 . 76 N 2 , and the stoichiometric chemical 
reaction for the system can be written as 

2H2 + O2 + 3.76N2 — »■ 2H2O + 3.76N2 

Calculations have been carried out for a grid with 
197 points in the circumferential direction and 152 
points in the normal direction. This grid was chosen 
based on the earlier work by Ahuja et al. [20] where 
the flow field was shown to be adequately resolved with 
this grid. For the present stoichiometric hydrogen-air 
mixture, the C-J velocity is Mach 5.11. Thus, for the 
Mach 6.46 case, the projectile speed is significantly 
above the detonation velocity of the mixture. If the 
free-stream velocity of the projectile is around the C- 
J detonation velocity of the mixture, unsteady flow 
phenomenon can occur. As given in Table 1, we shall 
be discussing five cases, two cases with Mach 5.11 and 
projectile diameters of 15 mm and 2.5 mm, and the 
remaining three cases with a Mach number of 6.46 and 
projectile diameters of 15 mm, 150 mm, and 225 mm. 
In all the calculations the residuals dropped by three 
orders in 12,000 iterations and then remained constant. 
Case (i): Mach 5.11 and projectile diameter of 15 
mm 

First the projectile diameter was chosen as 15 mm 
for the free-stream Mach number of 5.11. Figure 3 
shows the contour plot of temperature, and Fig. 4 
shows the contour plot of water mass fraction. The 
bow shock is very smooth, but the reaction front shows 
oscillations, which is more clearly seen in the inset 
enlarged view. The bow shock and the reaction front 
are separated from each other by the induction distance 
as is clearly seen in Fig. 5 which is the contour plot of 
density. The separation is minimum at the stagnation 
line and increases away from it because of low post- 
shock temperature away from the stagnation region. A 
comparison with Fig. 1 shows that all flow features are 
very well captured. The calculated frequencies of these 
oscillations [20] were found to be in good agreement 
with the experimentally observed frequency [5, 21]. 

Case (ii): Mach 5.11 and projectile diameter of 2.5 
mm 
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Next the diameter of the projectile was reduced 
to 2.5 mm while keeping the same ffee-stream Mach 
number of 5. 1 1. Other free-stream conditions were also 
kept the same. Figure 6 shows the temperature con- 
tours, and Fig. 7 shows the density contours. Both the 
shock and the reaction front are very smooth. This is 
much more clear from the water mass fraction contours 
shown in Fig. 8 with an inset enlarged view showing 
smooth reaction front. Thus, reducing the projectile 
diameter caused the instabilities to disappear. 

Case (iii): Mach 6.46 and projectile diameter of 15 
mm 

In this case the incoming Mach number was in- 
creased to 6.46, making it a superdetonative case, while 
keeping the projectile diameter of 15 mm as in case 
(i). Figure 9 shows the contour plot of density for 
Mach 6.46. The bow shock and the reaction front can 
be clearly seen in the figure. They are coupled with 
each other near the stagnation line and up to about 60 
degrees from the nose, at which point they start de- 
coupling from each other by the induction distance. 
This occurs because bow shock is almost normal near 
the stagnation line and the post-shock temperature is 
maximum. For Mach 6.46, a very small induction dis- 
tance occurs as a result of the post-shock temperature 
remaining significantly high up to some distance near 
the stagnation zone. Away from the stagnation line, 
the induction distance is increased as a result of de- 
creasing shock strength and post-shock temperature. A 
comparison with Fig. 2 shows that all the flow features 
are very well captured. Figure 10 shows the contour 
plot of temperature. A comparison with Fig. 3 shows 
that when the Mach number is increased from a C-J 
Mach number of 5.1 1 to a superdetonative Mach num- 
ber of 6.46, while keeping the same projectile diameter, 
the instabilities in the reaction front disappear. Thus, 
for the Mach 6.46 case with projectile diameter of 15 
mm, both the bow shock and the reaction front have 
a smooth profile. 

Case (iv); Mach 6.46 and projectile diameter of 150 
mm 

Next the projectile diameter was increased to 150 
mm while keeping the same Mach number of 6.46. 
Figure 1 1 shows the enlarged view of contour plot of 
water mass fraction for Mach 6.46, but with projectile 
diameter of 150 mm. The figure shows the waviness 
of the reaction front, indicating that the reaction front 
is turning towards instability. 

Case (v): Mach 6.46 and projectile diameter of 225 
mm 

Next, the diameter of the projectile was increased 
to 225 mm while keeping the same free-stream Mach 
number of 6.46 as in the previous case. Figure 12 


shows the enlarged view of the water mass fraction 
contours, and the periodic instabilities of the reaction 
front are clearly evident. Thus, increasing the projectile 
diameter from 15 mm to 150 mm and then to 225 mm, 
while keeping the same superdetonative Mach number 
of 6.46, causes the transition of a stable reaction front 
to an unstable one. The results of the above five cases 
have been summarized in Table 2. 


Table 2 

Summary of five test cases 


Moc 

Projectile 

Diameter 

D 

Profile of 
Reaction Front 

5.11 

(i) 15 mm 

(ii) 2.5 mm 

Unstable 

Stable 

6.46 

(iii) 15 mm 

(iv) 150 mm 

(v) 225 mm 

Stable 

Slightly 

Unstable 

Unstable 


In order to have a clear understanding of the origin 
and propagation of the instabilities, an x-t diagram 
for water mass fraction on the stagnation streamline 
between the bow shock and the reaction front is drawn. 

Figure 13 shows the x-t plot of water mass frac- 
tion along the stagnation line for the Mach 5.11 for a 
15 mm projectile diameter with an overlay of pressure 
to show the location of shock front. Fig. 14 shows 
the x-t plot of water mass fraction along the stagna- 
tion line for Mach 5.11, but with a projectile diame- 
ter of 2.5 mm. In the former case the reaction front 
clearly shows periodic oscillations, whereas the later 
case shows a smooth reaction front. The instabilities 
for Mach numbers lower than the C-J Mach number 
are due to the ignition delay. What causes these in- 
stabilities to disappear for the same Mach number of 
5.11 but lower projectile diameter will be explained 
in subsequent paragraphs. A similar trend which was 
observed for the superdetonative case of Mach 6.46, 
where increasing the diameter causes the instabilities to 
reappear, shall also be discussed. A qualitative study 
of the frequencies, shock structure, pressure, and stag- 
nation temperature associated with the projectile veloc- 
ities, shows that the reciprocal frequencies or periods 
of oscillation are equal to the induction time for the 
ignition of hydrogen-air mixtures. 

Many features of the non-steady flow regimes ap- 
pear to originate from near the stagnation region behind 
the normal segment of the bow shock. Also, there is 
clearly a relationship between the periodicity of these 
flow features and the magnitude of the chemical induc- 
tion time. The appearance and disappearance of these 
instabilities with various nose diameters have been ex- 
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plained with a wave interaction model as proposed by 
McVey and Toong [9]. Figure 15 shows a schematic 
x-t diagram for the wave interaction model for Mach 
5.11 and with projectile diameters of 15 mm and 2.5 
mm. For the diameter of 15 mm, the beginning of 
the cycle is shown at a time when the overtaking of 
the bow shock by a compression wave in the induc- 
tion zone takes place. Since chemical induction time 
is exponentially dependent on the gas temperature, the 
contact discontinuity (DCi) generated when the inci- 
dent compression wave overtakes the bow shock has 
an important effect on the location of a reaction front. 
Exothermic reaction takes place on the upstream side of 
the contact discontinuity (DCi) before beginning on the 
downstream side due to the difference in the gas tem- 
perature on the upstream and downstream sides of the 
contact discontinuity. It is this beginning of exothermic 
reaction on only one side of the contact discontinuity 
which constitutes the formation of a new reaction front. 
The ratio of chemical induction time across the discon- 
tinuity is a measure of where the new reaction front will 
be located relative to the original reaction front. Since 
the gases are hot on the upstream side of the contact 
discontinuity, they begin to react, generating compres- 
sion waves which propagate upstream and downstream. 
This is shown in the figure by (CW 2 ) and (CW 3 ). At a 
somewhat later time, the contact discontinuity reaches 
the position of the original reaction front, extinguish- 
ing the reaction at this front and generating rarefac- 
tion waves (RWj) and (RW 2 ). The upstream-facing 
compression wave (CW 2 ) eventually overtakes the bow 
shock and thus results in the generation of another con- 
tact discontinuity (DC 2 ). Presence of both rarefaction 
and compression waves in the induction zone makes it 
possible for a periodic wave interaction process. 

It is clear from the above discussion that a com- 
pression wave moving towards the bow shock is over- 
taken by a rarefaction wave in each cycle of the wave- 
interaction process. During this process the compres- 
sion wave gets attenuated by a rarefaction wave. If 
Msi is the Mach number of the incident compression 
wave (CWi) or weak shock, then (Msij^ax is generally 
of the order of 1.01. A wave interaction process cannot 
be initiated by an unattenuated compression wave since 
the Mach number ( Msi) of such a wave is much greater 
than (Msi)max* On the other hand, if a sufficiently weak 
compression wave does initiate the wave interaction 
process, the attenuation of the first and each subse- 
quent compression wave by a rarefaction wave must 
be sustained in order that the Mach number of each 
reaction shock, which finally overtakes the bow shock, 
be less than (Msi)n^. When this condition is satisfied, 
then only the cyclic process is feasible. The validation 
of this phenomena with a shock-fitting method shall be 
presented in future studies. 

The importance of a compression wave reflecting 
off the projectile nose has not been clarified in the 


McVey and Toong model. The compression wave 
which travels towards the projectile reflects from it and 
finally overtakes the bow shock. This overtaking of 
the bow shock by the reflected compression wave from 
the projectile takes place almost at the same time that a 
newly formed compression wave from the new reaction 
front has overtaken the bow shock. Consequently, 
this contributes to strengthening the compression wave 
coming from the reaction front. 

Thus, when we compare the unstable Mach 5.11 
case with a projectile diameter of 15 mm with the sta- 
ble case of a projectile diameter of 2.5 mm, it is clear 
that the 15 mm diameter case has greater shock stand- 
off distance than that with a diameter of 2.5 mm. Thus, 
the reflected compression wave from the projectile nose 
coming back to overtake the bow shock to strengthen 
the compression wave originating at the reaction front, 
would have become quite weak. This is because it un- 
dergoes a large number of interactions with other com- 
pression waves between the reaction front and the pro- 
jectile nose. As shown in the figure, the strength of the 
shock P5/P3 is less than the strength of the original re- 
flected shock which is P 2 /P 1 . Now, if the shock stand- 
off distance is larger, the reflected compression wave 
would encounter a greater number of interactions, and 
with every interaction, its strength would decrease. Al- 
though this reflected compression wave strengthens the 
compression wave originating at the reaction front, the 
Mg is still less than (Msi)max and, therefore, the periodic 
instabilities are sustained. On the other hand, when the 
projectile diameter is 2.5 mm, the reflected compres- 
sion wave has to travel a very short distance before it 
reaches the reaction front to strengthen the compres- 
sion wave originating there. Moreover, because of the 
short shock stand-off distance, it does not encounter 
many interactions with other compression waves. Al- 
though P3/P2 of the reflected wave is less than P2/P1 of 
the compression wave reaching the projectile, it is still 
very strong by the time it reaches the reaction front. 
Thus, when it overtakes the bow shock, it strengthens 
the compression wave much more than the preceeding 
15 mm diameter case. This causes Ms of the compres- 
sion wave to be much greater than (Msi)max, and this 
leads to unattenuated compression waves and, hence, 
periodicity disappears. 

Again, when referring to the Mach 6.46 case with 
a projectile diameter of 15 mm, a stable reaction front is 
observed, whereas increasing the diameter to 225 mm 
causes the instabilities to reappear. For the 225 mm 
case where the shock stand-off distance is much higher 
than the 15 mm case, the shock wave reflected off the 
projectile body becomes much weaker before strength- 
ening the compression wave generated at the new re- 
action front. Thus, the compression wave generated at 
the new reaction front has Mj less than (Mji )max and, 
therefore, a cyclic process is sustained. On the other 
hand, for the projectile diameter of 15 mm where the 
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shock stand-off distance is much smaller, the reflected 
compression wave from the projectile body strengthens 
the compression wave generated at the new reaction 
front much more and, therefore, leads to unattenuated 
compression waves. Consequently, the cyclic process 
for the 15 mm diameter and Mach 6.46 case cannot 


be sustained. 

CONCLUSIONS 

Numerical studies with different projectile diame- 
ters and Mach numbers indicate that Chapman-Jougei; 
velocity is not the only parameter for triggering the 
instabilities. The other parameters are induction time, 
reaction rate constant, projectile diameter, activation 
energy and heat release. Ihe reason fcx the appearance 
and disappearance of instabilities could be attributed to 
attenuated or unattenuated compression waves which 
depend iipon the strength of the reflected compression 
wave from projectile body, the strength of which in 
turn depends upon the shock stand-off distance and, 
hence, the projectile diameter. Thus, it is concluded 
that a subdetonative case can also lead to a stable re- 
action front by having an appropriate small diameter 
projectile, and a superdetonative case can also lead to 
an unstable reaction front by having a larger diameter 
projectile. 
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Figure 1 Shadowgraph of a spherical nose projectile moving at 
Mach 5.11 into a premixed stoichiometric hydrogen-air mixture. 



Figure 2 Shadowgraph of a spherical nose projectile moving at 
Mach 6.46 into a premixed stoichiometric hydrogen-air mixture. 
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Figure 3 Contour plot of temperature for Mach 5.11 and projectile diameter 15 mm. 
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Figure 4 Contour plot of water mass fraction for Mach 5.11 and projectile diameter 15 mm. 
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Figure 5 Contour plot of density for Mach 5.11 and projectile diameter 15 mm. 
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Figure 6 Contour plot of temperature for Mach 5.11 and projectile diameter 2.5 mm. 
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Figure 7 Contour plot of density for Mach 5.11 and projectile diameter 2.5 mm. 



Figure 8 Contour plot of water mass fraction for Mach 5.11 and projectile diameter 2.5 mm. 
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Figure 9 Contour plot of density for Mach 6.46 and projectile diameter 15 
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Figure 10 Contour plot of temperature for Mach 6.46 and projectile diameter 


mm. 


15 mm. 
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Figure 11 Enlarged view of contour plot of water mass 
fraction for Mach 6.46 and projectile diameter 150 mm. 
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Figure 12 Enlarged view of contour plot of water mass 
fraction for Mach 6.46 and projectile diameter 225 mm. 
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Figure 13 x-t Plot of water mass fraction along stagnation 
streamline for Mach 5.11 and Projectile diameter 15mm 
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Figure 14 x-t Plot of water mass fraction along stagnation 
streamline for Mach 5.11 and Projectile diameter 2.5 mm 
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